Exciton-phonon effects in carbon nanotube optical absorption 
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We find that the optical properties of carbon nanotubes reflect remarkably strong effects of exciton- 
phonon coupling. Tight-binding calculations show that a significant fraction of the spectral weight 
of the absorption peak is transferred to a distinct exciton+phonon sideband, which is peaked at 
around 200 meV above the main absorption peak. This sideband provides a distinctive signature 
of the excitonic character of the optical transition. The exciton-phonon coupling is reflected in a 
dynamical structural distortion, which contributes a binding energy of up to f 00 meV. The distortion 
is surprisingly long-ranged, and is strongly dependent on chirality. 
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The optical properties of carbon nanotubes are cur- 
rently the focus of intense experimental and theoretical 

attention HEHEEH0,lEHEl[IlE2, and even 
single nanotube electro-optical devices have been demon- 
strated [Isl H3. 0|. Most of the experimental results 
have been discussed and analysed in terms of interband 
transitions. In particular, nanotube bandgaps were de- 
termined in this way. However, theoretical calculations of 
the optical spectra 0, 0, 0] suggest that the observed 
transitions correspond to exciton energies, not interband 
transitions. Moreover, the exciton binding energies are 
anomalously large in nanotubes 0, El El U3 cor ~ 
responding to a substantial fraction of the bandgap, so 
the optical transition energy is quite different than the 
bandgap. 

If correct, these theoretical results require a reevalua- 
tion of our current picture of nanotube electronic struc- 
ture. Device properties are particularly sensitive to the 
bandgap, because of the central role of tunneling in most 
nanotube transistors [l^. Thus a definitive test of the 
excitonic interpretation of optical transitions is needed. 
Fine structure in the optical data can provide a dis- 
tinctive signature, facilitating comparison between differ- 
ent models to unambiguously verify the role of excitonic 
transitions. Recent photoconductivity excitation spec- 
tra [rH Hi| show sidebands at about 200 meV above to 
the main absorption peaks, suggesting the involvement 
of phonons. Fluorescence excitation spectroscopy also 
shows structure at the phonon energy, allo wing more 
detailed analysis with finer energy resolution 19] . But 
while the interaction of single electrons with phonons 
has been studied extensively, we know of no studies of 
exciton-phonon interactions in nanotubes. 

We therefore investigate theoretically the role of 
electron-phonon coupling in the optical spectra of nan- 
otubes, including a comparison of between excitons and 
free-carrier transitions. We find surprisingly strong 
phonon effects in the excitonic spectra. Dynamical ef- 
fects lead to the transfer of a significant fraction of the 
spectral weight from the exciton absorption peak (zero 
phonon line) to a phonon sideband peaked at around 



200 meV (Fig. QJi and Fig. Hk-b), consistent with the 
experimental observations. In contrast, our calculations 
for electron-phonon interaction without excitonic bind- 
ing (Fig. Eh) do not show a distinct phonon sideband. 
These results clearly suggest that the optical data must 
be interpreted in terms of excitonic transitions, and do 
not provide a direct measure of the bandgap. 

We find that the fraction of the intensity transferred 
to the phonon sideband is inversely proportional to the 
nanotube diameter. The exciton-phonon binding energy 
is also unexpectedly high, about 60-100 meV. The asso- 
ciated lattice distortions exhibit an intriguing structure, 
they extend far beyond the exciton itself [llj and reverse 
sign for different nanotube chiralities (Fig. ^Jd-c). The 
effect of the dielectric environment on the spectra is also 
discussed. (We treat individual nanotubes, as in fluores- 
cence excitation spectroscopy. Nanotube bundles are less 
suitable for measuring detailed lineshapes.) 

In the absence of exciton-phonon coupling, emission 
or absorption of a photon involves an exciton of total 
wavevector q$ corresponding to the photon momentum, 
hereafter approximated as = 0. The exciton-phonon 
coupling mixes this exciton with phonons and with ex- 
citons of other q, such that the total exciton+phonon 
momentum is conserved. The finite-g exciton wavefunc- 
tion can be found from the solution of the Bethe-Salpeter 
Equation (BSE) [20( in the basis of a tight-binding Hamil- 
tonian [21| , analogous to the q = case 01 : 

l^) = E^ u U^i GS >- (!) 

k 

Here Af q is the eigenvector of the <S"s state of BSE so- 
lution; u\ +q (vk) creation (annihilation) of an electron 
in the conduction (valence) band acting on the ground 
state |GS) = Ylk v k l vac )- The indices k and q each label 
both the continuous ID wavevector along the tube axis 
and the discrete circumferential wavevector. 

We model the electron-phonon interaction by the Su- 
Schrieffer-Heeger (SSH) model 113], with matrix element 
t = to — gSu dependent on the change of the nearest 
neighbor C-C distance (Su), where to = 3 eV. We take the 
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FIG. 1: (color online), (a) Absorption spectrum Eq. |gj 
for (17,0) tube and e = 2 (solid red curve and red delta- 
function). Dotted line shows electronic contribution to side- 
band energy (with 10 meV Gaussian broadening). Phonon 
contributions are quite narrow and are shown schematically 
as delta-functions (height ~ J 1 / 2 ). (b) Long-range tail of the 
atomic displacements with respect to the hole position (at ori- 
gin) for a (16,0) tube: radial displacements (labelled 5R, in 
blue), and displacements parallel to the axis (labelled 5z, in 
red), (c) Same for (17,0) tube, (d) Short-range azimuthal and 
axial distortions, shown on "unwrapped" carbon nanotube by 
displaying the displaced atom positions (red circles) together 
with undisplaced positions (white circles) , with displacements 
magnified by a factor of 230 for visibility. Hole position is solid 
circle at center. 
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FIG. 2: (color online). Absorption spectra in (17,0) tube 
calculated with and without electron-phonon couplings (solid 
red and dashed black curves respectively) for (a) e — 2, (b) 
e = 4, (c) e — > oo. The zero of energy here corresponds to the 
onset of the first optically active exciton, which has zero width 
and is shown by the vertical arrows. The onset of continuum 
(bandgap) transitions is at 0.48 eV, 0.19 eV, and in (a), 
(b), and (c) respectively. The width of the second exciton is 
finite due to: (1) coupling with the first band electron-hole 
continuum to give lifetimes of t cc = 10 fs and 15 fs for e =2 and 
4 respectively; and (2) coupling with phonons to give r p h=90 
fs and 33 fs for e =2 and 4 respectively. Delta functions 
in Eq. JSJ were Gaussian broadened with width of 20 meV, 
except in (c) we use 6.5 meV to emphasize that the width of 
the continuum peak is intrinsic. 



electron-phonon coupling constant to be g = 5.3 eV/A as 
predicted theoretically for a related molecular problem 
[23j, consistent with fits to the Peierls gap in conjugated 
polymers [24j. 

After Fourier transformation the intraband SSH 
Hamiltonian has the form: 



Hd-ph = M k q ( v l+ q v k - u\ +q u k ){a qil 

kqft 



+ al qil ), (2) 



_1 / 2 is momentum dependent electron- 
,t 



where M£ oc gN 

phonon coupling; a|_ gM is a phonon creation operator 
with wavevector — q and phonon band index fi = 1...6; 
and N is the number of primitive unit cells, each con- 
taining two carbons. For the phonon spectrum we used 
a force-constant model similar to Saito et al. |25j . 

The electron-phonon Hamiltonian mixes the optically 
active q=0 exciton with finite-g excitons in combination 
with phonons of wavevector —q: 



S'q/J, 



qH 



S' 



= E M "X*«o + ^ +9 ,o) (3) 
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Here the orthogonality relation of the BSE solution 
(Es%^f'g = <W) nas been used to derive exciton- 
phonon coupling amplitudes £?^f . 

The wavefunction and the spectral line shape can be 
evaluated in second order perturbation theory for the 
lowest optically active exciton s: 
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The calculated absorption spectrum Eq. (@J is shown 
in Fig. ^ (solid red curve) for a (17,0) tube and e = 
2, where e is the dielectric constant of the embedding 
medium Because of exciton-phonon coupling, the 

main absorption peak looses 8% of it's spectral weight 
to the sideband, which corresponds to the continuum of 
finite-q excitons plus phonon of wavevector —q. Most 
of the transferred spectral weight goes to the prominent 
sideband at about 210 meV above the zero phonon line, 
with 4% of the spectral weight falling between 100 and 
300 meV. The SSH Hamiltonian has little coupling to the 
low frequency modes and we find a much weaker replica 
at the radial breathing mode (RBM) frequency [i.e. 20 
meV for a (17,0) tube]. 

To understand the spectrum in more detail, we decom- 
pose the sideband energy into contributions from phonon 
energy and exciton dispersion. Specifically, for each tran- 
sition within 300 meV above the zero-phonon line, we 
project out the energy contribution from phonons vs elec- 
tronic excitation (due to the admixture of finite-g exci- 
tons). The phonon contribution to the sideband exhibits 
three peaks Fig. corresponding to the longitudinal 
optical (LO) phonon band edges at the K and V points 
of the graphene Brillouin zone, and the radial breath- 
ing mode. The K phonon dominates. It has stronger 
coupling; and more importantly, it mixes exciton bands 
[2fij |. allowing absorption by the "dark" (dipole forbid- 
den) bands [ll|,|27j. The energy difference between the 
sideband peak position and optical phonon frequency is 
due to the exciton "recoil energy" , i.e. the energy from 
finite exciton q, which contributes ~30 meV here. 

In the presence of an exciton the nanotube distorts dy- 
namically. While the electron remains within roughly 2 
nm of the hole for the cases shown in Fig. ^ the struc- 
tural distortions are far more long-ranged. The distor- 
tions calculated using the wavefunction Eq. Q are shown 
in Fig.^ 3 andQJ; for (16,0) and (17,0) tubes respectively, 
plotting the atomic displacements relative to the posi- 
tion of the hole. The breathing distortions decay expo- 
nentially away from the exciton, with a decay length of 
A5 = 13 nm. We find that is proportional to the tube 



diameter. The sign of the breathing distortion depends 
on chirality indices (n,m): positive for mod(n-m,3)=l 
and negative for mod(n-m,3)=2. From the sign of the 
exciton-phonon matrix elements, we expect a reversal of 
the signs of the breathing distortions for the second exci- 
ton. The displacements parallel to the tube axis decay far 
more slowly even than this, oscillating with a wavevector 
A z ~ 80 nm, which appears to be insensitive to the tube 
diameter. The short-range distortions near the exciton 
are shown in Fig. ^i. 

The spectra for higher-energy excitons (Eq > Eq) can- 
not be obtained from Eq. J3J, because the denominator 
Eg + huj- qt i — Eq can be arbitrarily close to zero and the 
perturbation theory breaks down. Toyozawa showed |2sj 
that an exact solution for the absorption spectra has a 
form similar to the perturbation theory expression, with 
an energy dependent lifetime and polaronic shift in the 
energy denominator. We approximate Toyozawa's solu- 
tion by evaluating the lifetime broadening in the Random 
Phase Approximation (RPA): 



I(w) 



E fs 

S 



(hcj-E^ + Tsicuy 



TsH - 7rE|^'| 6(Huj-E( -Huj^), (5) 

S'qfJ, 

where fs is the oscillator strength of the S exciton [llj . 
We checked that the Toyozawa solution obtained with a 
self consistent Born approximation neglecting the k de- 
pendence of the self energy Eg (E) = As(E)+Ts(E) does 
not change the RPA result Eq. (0. If huj — Eq 3 in Eq. 
(0 then Ts(Eq)/H equals half of the reciprocal lifetime 
of S'-exciton due to the scattering by phonons. The bind- 
ing energy shift for the 5-exciton in RPA approximation 
is: 



S'qn 



V{E» - Ef ~ 



(6) 



where V denotes the principle part. 

The absorption spectra calculated with and without 
exciton-phonon coupling are shown in Fig. for the tube 
embedded in dielectric e = 2, 4, and for the free electron- 
hole pair absorption (equivalent to the limit e — > oo). In 
the absence of exciton-phonon interactions there are two 
strong exciton absorption peaks, each followed b y th e cor- 
responding continuum of intraband absorption |ll| . The 
first exciton has zero width in this approximation, and 
is shown by a vertical arrow in Fig. and Eb- Both 
exciton lines show a distinct phonon sideband about 200 
meV above the main absorption line. (For e = 4 the exci- 
ton binding energy is about 0.2 eV, so the band-to-band 
absorption (dotted line) also contributes to the sideband 
intensity in Fig. |2b-) In contrast, there is no distinct 
phonon peak associated with band-to-band absorption 
(Fig. |2t). Thus the phonon peak provides a clear sig- 
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FIG. 3: (color online), (a) The phonon binding energy, and 
(b) fraction of the spectral weight transferred to the exci- 
ton+phonon sideband, for the first exciton in zig-zag tubes, 
versus tube diameter for e — 2 (black circles) and for e — 4 
(blue squares) respectively. The curves (black solid for e = 2 
and blue dashed for e = 4) are fits Atj + B b jjd, where 
A b = 47 (meV), B b = 30 (meVnm) and Ai = 1.7 (%), 
Bi = 10 (%nm) for e = 2. 



nature of whether or not the absorption is excitonic in 
nature. 

The second exciton can decay into free electron-hole 
pairs of the first band via Coulomb interaction or by emit- 
ting a phonon, giving rise to a finite lifetime of the second 
exciton resonance. The electronic and phonon contribu- 
tions are given in the figure caption for a (17,0) tube for 
two values of e. The actual value of the lifetime, as well 
as other details of the spectrum, are very sensitive to the 
lineup of the resonance peak position with the onset of 
the first band continuum, and thus are also sensitive to 
the tube radius and dielectric environment. Therefore 
experimental measurements of the first exciton are easier 
to compare with theory. 

The total exciton-phonon binding energy calculated 
from Eq. © has the largest contribution from the cou- 
pling to the higher energy states, which have smaller op- 
tical spectral weight due to the energy denominator in 
Eq. ©. The dependence of the binding energy on the 
diameter d is shown in Fig. Ol, along with a phenomeno- 
logical fit, for the first exciton in different zig-zag tubes, 
for e = 2 and 4. The second exciton binding energy is 
always larger than the first exciton, by 15-30%, due to 
the larger effective mass of the former. The fraction of 
the spectral weight transfer for the first exciton Eq. (0J 
is shown on Fig. |3J>, along with a fit for the same values 
of e. 

In conclusion, we calculate exciton-phonon effects in 
the absorption spectra of carbon nanotubes, predicting a 
significant spectral weight transfer from the main (exci- 
tonic) absorption peak to a phonon sideband at around 
200 meV above the zero phonon line, for both the first 
and second excitons. Comparable structure is seen in 



recent experiments 0,0,0]. In contrast, for band-to- 
band transitions we find that there is no distinct side- 
band; instead the phonon structure represents only a 
slight perturbation of the spectrum. We therefore be- 
lieve that the phonon sideband provides direct experi- 
mental evidence, that the optical transition is excitonic 
in nature and cannot be used as a direct measure of the 
bandgap. 
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